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The observed flux of ultra-high energy (UHE) cosmic rays (CRs) guarantees the presence of high- 
energy cosmogenic neutrinos that are produced via photo-hadronic interactions of CRs propagating 
through intergalactic space. This flux of neutrinos doesn't share the many uncertainties associated 
with the environment of the yet unknown CR sources. Cosmogenic neutrinos have nevertheless 
a strong model dependence associated with the chemical composition, source distribution or evo- 
lution and maximal injection energy of UHE CRs. We discuss a lower limit on the cosmogenic 
neutrino spectrum which depends on the observed UHE CR spectrum and composition and relates 
directly to experimentally observable and model-independent quantities. We show explicit limits 
for conservative assumptions about the source evolution. 

PACS numbers: 98.70.Sa,95.55.Vj 



I. INTRODUCTION 

Cosmogenic neutrinos are produced when UHE CRs 
interact with the cosmic radiation background while 
propagating between their sources and Earth. The 
frequent interactions with the cosmic microwave back- 
ground (CMB) limits the propagation of nucleons with 
energies greater than Eqzk — 40 EeV to within a few 
100 Mpc and is responsible for the so-called Greisen- 
Zatsepin-Kuzmin (GZK) cutoff of extra-galactic pro- 
tons [312]. Mesons produced in these interactions quickly 
decay and produce an observable flux of cosmogenic (or 
GZK) neutrinos 3J. In fact, the observed spectrum of 
CRs extending up to energies of a few 100 EeV shows a 
suppression above ~ -Egzk with high statistical signifi- 
cance 4, 5 . This could be an indication that protons are 
dominating the flux at these energies. In this case the 
flux of cosmogenic neutrinos is typically large. 

However, the experimental situation is less clear. Mea- 
surements of the elongation rate distribution of UHE CR 
showers indicate a transition of their arrival composi- 
tion from light to heavy within 4-40 EeV O [7]. If a 
heavy component dominates also at higher energies the 
prospect for cosmogenic neutrino production is "disap- 
pointing" [8] or at least less favorable than for the pro- 
ton scenario. A crucial uncertainty of this scenario is 
the maximal injection energy of the nucleus with mass 
number A; as long as E max ^S> AEqzk, even this sce- 
nario will produce an appreciable amount of cosmogenic 
neutrinos [5]. If this condition is not met interactions 
with the subdominant cosmic photon background from 
the optical/infra-red will still contribute to the cosmo- 
genic neutrino flux. We will use the estimate of Ref. [10] 
for our calculation. 

The IceCube neutrino observatory has reached the sen- 
sitivity for the detection of optimistic cosmogenic neu- 
trino fluxes [11] . In the case of a non-observation it is 
of interest to know a lower limit on the various source 
emission possibilities for their definite exclusion. Lower 



cosmogenic neutrino flux limits have already been dis- 
cussed in the context of proton-dominated scenarios via 
a deconvolution of early Auger data [T2] • We will discuss 
in this article updates of these lower limits and exten- 
sions to more general assumptions for the source distri- 
bution and chemical composition. Similar to Ref. |12j we 
will not attempt to construct a specific source emission 
model that fits the Auger spectrum and elongation rate 
distribution but we will derive the limits directly from the 
observed composition measurement and spectrum. From 
this we can derive a strict lower limit on the cosmogenic 
flux. 



II. COSMIC RAY PROPAGATION 

The propagation of UHE CR nuclei is affected 
by photo-disintegration [13l [14], photo-hadronic 
interactions [TS], Bethe-Heitler pair production [TB] and 
red-shift losses due to the expansion of the Universe. It 
is convenient to consider a homogenous and isotropic 
distribution of CR sources and derive the observed CR 
from the co-moving number density Y± = rij/(l + z) 3 as 
a solution to a set of Boltzmann equations [T7] , 

Y i = d E {HEY i ) + d E (b i Y i ) 

-Tl ot Y t + J2 J d^7^+A- (1) 
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The cosmic expansion rate H(z) follows the usual "con- 
cordance model" dominated by a cosmological constant 
with f^A ~ 0.73 and a (cold) matter component, £! m ~ 
0.27 with H 2 (z) = H$ [O ro (l + z) 3 + fi A ], normalized to 
its value today of H ~ 72 kms" 1 Mpc" 1 [H]. The first 
and second terms on the r.h.s. of Eq. (JlJ describe, respec- 
tively, red-shift and other continuous energy losses (CEL) 
with rate b = —dE/dt. In the following we will treat 
Bethe-Heitler pair production as a CEL process [16] . The 
third and fourth terms describe more general interactions 
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involving particle losses (i — > anything) with total inter- 
action rate r* ot , and particle generation of the form j — >• i 
with differential interaction rate 7.^ . The last term on the 
r.h.s., Ci, corresponds to the emission rate density of CRs 
of type i per co-moving volume. The detailed description 
of the interaction rates and their scaling with red-shift 
has been discussed in our previous publications [T7[ Q]|] . 

We first discuss the case of proton sources. The flux of 
cosmogenic neutrinos today [z = 0) depends on the co- 
moving number density of protons at all red-shifts and 
can be approximated as [17] 
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x / d£ p j pI/ (z',£ p ,(l + z')E u )Y p (z',£ p ), (2) 

where £ p is the solution to the differential equation £ p — 
—H£ p — &bh(z, £ p ) with initial condition £ p (0, E p ) = E p . 
The co-moving number density of protons can be written 
as 
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where the effective source term is defined as 



£p,e«{ z , E p ) 

= C p (z,E p ) + I d£ p 'y pp (z,£ p ,E p )Y p (z,£ p ) . (4) 



III. MINIMAL NEUTRINOS FROM PROTONS 



A minimal contribution to the flux of cosmogenic neu- 
trinos can be estimated as follows. As a first step we 
approximate the UHE CR spectrum measured by Auger 
via the phenomenological fit given in Ref. [20]. This fit 
is shown in Fig. [T] as a dashed-dotted line together with 
recent data of Auger, HiRes [J] and the Telescope Ar- 
ray [21] (TA). Note, that the normalization of the Auger 
data is lower by a about a factor two than HiRes and TA 
and hence cosmogenic neutrinos derived from this data 
are the lowest. 

Whereas the spectrum of UHE CRs is dominated by 
closeby sources, the neutrino flux receives contributions 
up to the Hubble scale. The overall flux will hence in- 
crease for an increasing number of sources with red-shift. 
We assume that redshift evolution decouples from the 



source emission spectrum, i.e. C p (z,E) — H(z)Q p (E) 
and we consider two scenarios for the source evolution 
T-L(z). In the most conservative case we assume source 
contributions within redshift z max = 2 with no source 
evolution, i.e. Ho = Q(z max — z). A more realistic sce- 
nario assumes a source evolution following the star for- 
mation rate. We will use the estimate [221 1231 
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(l + z) 3A 
Ni (1 + z)-°- 3 
Ni N 4 (1 + z)- 3So 



z < 1 , 

1 < z < 4, 

z > 4, 



(5) 



with normalization factors, Ni = 2 3 7 and N4 — 5 3 2 . 
Since we assume conservative choices of the source evo- 
lution the associated cosmogenic neutrino flux can be re- 
garded as lower limits on the expected cosmogenic neu- 
trino flux. 

In the following we will derive approximate solutions 
to Eqs. ([2]) and Q using an iterative scheme. For the 

iteration start we choose Q p °\e p ) — (Hq + Oeoq + 
r )4TrJ^(E p ), where b and T are the energy loss and 
interaction rate, respectively, at redshift z = 0. The it- 
eration step is then given by 



Q p n+1 \E P ) = 4nJ$g(E p )/r,W(E p ) , 



(6) 



where we use the phenomenological fit of Ref. [20] for 
Jcr(E) and introduce the effective survival distance 
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We continue this iteration until the relative correction 
Ei(Q%? l) /Q%i ~ !) 2 st °P s t0 decrease or a maximal 
(sufficiently large) iteration step is achieved. This com- 
pensates for numerical instabilities. 

In Figure [l] we show the resulting cosmogenic neutrino 
flux for this procedure for the two evolution scenarios. 
The limit for the SFR evolution agrees well with that de- 
rived from a deconvolution analysis in Ref. [12] . We also 
indicate in this plot the sensitivity of IceCube [11] and 
the proposed Askaryan Radio Array (ARA) [23]. Three 
years of observation with the 37 station configuration of 
ARA ( ' ARA-37" ) is sufficient to reach the proton emis- 
sion model for the SFR case. In the case of no source 
evolution this scenario is reached after ten years. As we 
already emphasized, this result depends on the absolute 
normalization and/or energy calibration of the observed 
UHE CR spectrum. For a normalization to HiRes and 
TA data we expect our limits to scale up by about a 
factor 2. 
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FIG. 1: Minimal flux of cosmogenic neutrinos assuming dom- 
inance of protons above 4 EeV. We show the results without 
source evolution (dotted) and assuming source evolution ac- 
cording to the star formation rate (solid). Also shown are the 
projected sensitivities of IceCube (10 years) and the ARA-37 
(3 years) as dashed lines. The thick dashed-dotted line shows 
the approximation of the Auger spectrum above the ankle. 
For comparison, we also show the bestfit cosmogenic neutrino 
flux (green solid line) from Ref. [24] (£ min = 10 18 5 eV) in- 
cluding the 99% C.L. (green shaded area) obtained by a fit to 
the HiRes spectrum. 



IV. GENERALIZATION TO HEAVY NUCLEI 

The case of a more general scenario including UHE 
CR sources of heavy nuclei is more complicated. The 
chemical composition observed at Earth is the result of 
rapid photo-disintegration in the radiation background 
and there is no simple connection to the source com- 
position. However, since photo-disintegration conserves 
the energy per nucleon we can derive a lower neutrino 
limit by tracking the leading (heaviest) nucleus back to 
its source starting from a composition A a and Z a inferred 
from UHE CR observations. 

The parent nuclei during this back-tracking are at least 
as heavy as the observed mass composition. For instance, 
a single helium nucleus in the observed spectrum might 
be produced via the production chain 10 B — > 9 Be(+p) — > 
4 He(+ 4 He + p) from the source. The parent nuclei in 
each step of this chain determine the interaction and en- 
ergy loss rates during propagation. For a lower limit on 
the cosmogenic neutrino flux we have to minimize the 
emission rate density of the UHE CR nuclei associated 
with their cascades in the CMB. This corresponds to a 
maximal survival probability of nucleons. Hence, we can 
derive a strict lower limit with the assumption that the 



back-tracking of the nuclei is indefinite, i.e. we assume 
no upper limit on the atomic mass number in the nuclei 
cascades. 

Photo-disintegration that drives the cascades competes 
with photo-hadronic interactions and Bethe-Heitler en- 
ergy loss. To first order, a photo-hadronic interaction of 
the nucleon with energy E, charge Z and mass number 
A can be approximated via the interaction rate of the 
free proton as Taj(E) ~ AT pl (E /A). Hence, the inter- 
action rate per nucleon of the parent nucleus is the same. 
Energy loss via Bethe-Heitler pair production, however, 
scales as bA-y{E) ~ Z 2 b pl (E /A) and the effective energy 
loss per nucleon scales as Z 2 /A. Again, for a maximal 
survival probability of the nucleons and hence a mini- 
mal emission rate density of the sources, we assume a 
minimal Bethe-Heitler energy loss of the nucleons. This 
corresponds to the energy loss of a nucleus with charge 
Z a and atomic mass number A a associated with the ob- 
served composition. 

In summary, a lower limit on the cosmogenic neutrino 
flux can hence be derived by the same Eqs. ^ and 
([3]) where we now replace the continuous energy loss by 
its minimal contribution b m - m (z,E) ~ (Z 2 / A )b^(z , E) , 
where £>bh correspond to the energy loss of a free pro- 
ton. The photohadronic interaction of the nucleons is 
given by the average interaction of protons and neutrons. 
The total number of nucleons per nucleon energy depends 
on the observed (or inferred) mass composition of UHE 
CRs. Assuming a single component we have the relation 
E N J N (E N ) = A E C rJcr(E C r) with E N = E C r/A q or 
Jn(En) — A J C r(Ecr)- 

In the left panel of Fig. [2] we show the minimal cos- 
mogenic neutrino fluxes for the case of helium, nitrogen, 
silicon and iron dominance of the Auger spectrum. The 
level of these fluxes is not in reach of present or future 
neutrino observatories. However, cosmogenic neutrino 
fluxes strongly depend on the maximal injection energy 
of the sources. We conservatively assume for our method 
that the maximal energy does not exceed the observed 
energy of UHE CRs. However, it is in principle possible 
that these models produce detectable fluxes of cosmo- 
genic neutrinos [9] if the maximal energy significantly 
exceeds A x Eqzk- We will briefly discuss this in the 
following section. 

We can also generalize our method to the case of a 
mixed compositions, which is indicated by the Auger CR 
elongation rate distribution. For instance, if /j(-E'cR) de- 
notes the fraction of nuclei with mass Aj at CR energies 
Ecr the mean mass number is given by 

J N (E N ) ~ AlfiiAiE^JcniAiE^) . (8) 

Hence the minimal cosmogenic neutrino flux in this 
case is J™' ,n (E u ) = Ei^O^). whcrc the individual 
jmm are derived m the same way as before but using 
fi(EcR) Jcr(Ecr) a s the input spectrum. As an example 
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FIG. 2: Minimal flux of cosmogenic neutrinos for a mixed composition. Left panel: Minimal flux of cosmogenic neutrinos 
assuming dominance of protons, helium, nitrogen, silicon or iron in UHE CRs above 4 EeV. We show the results without source 
evolution (dotted) and assuming source evolution according to the star formation rate (solid). Right panel: The contribution 
of protons (red lines) in a mixed composition scenario assuming 100% (upper line), 10% (middle line) and 1% (lower line) 
proton contribution (black lines) at 100 EeV. 



we show in the right panel of Fig.[2]the lower limit associ- 
ated with protons in a multi-component model, where we 
decrease the proton contribution at 100 EeV to 10% (a = 
1) and 1% (a = 2) using f p = 1 - (1+ (E/lO^eV)- 01 )- 1 
with / A = 1 - f p . 



V. OPTIMISTIC COSMOGENIC NEUTRINOS 



Predictions of the cosmogenic neutrino spectra are very 
sensitive to the maximal energy of UHE CR nuclei. In the 
following we will briefly discuss "optimistic" predictions 
that assume that the maximal energy of CR nucleons is 
much larger than the GZK cutoff, i.e. Ecr/A ^> Eqzk- 
For the discussion it is convenient to introduce the energy 
density (eV cm" 3 ) of the GZK neutrino background at 
redshift z defined as 



w GZ K = J &E V E V Y V {E V ). 



(9) 



From the Boltzmann equations we can derive the evo- 
lution of the energy density as 

uj G zk + Huj gzk = J2 [ dEb itGZK (z,E)Y i {z,E), (10) 

where bi^QZK(E) — 0.2ET 77r (E/Ai) is an approximation 
of the energy loss of the nuclei into GZK neutrinos [9] . 

The UHE CR interactions with background photons 
are rapid compared to cosmic time-scales. The energy 



threshold of these processes scale with redshift z as 
AiE t h/ '(1+z) where i? t h ^ Eqzk is the (effective) thresh- 
old today. We can therefore approximate the evolution 
of the energy density as 



wgzk + Hujqzk 



ZKirHjz) 
4(1 + K v 



J dEEQ z (E), (11) 



AiE th /(l+z) 



where is the ratio of charged to neutral pions pro- 
duced in j»7 interactions. Assuming a power-law emission 
rate density Qi(E) oc E~ 11 with sufficiently large cutoff 
E max ^> E t h we see that cosmic evolution enhances the 
GZK flux as 

3^ (AiE^QiiAiE^) 



(12) 



where the last term assumes ji > 2 and the effective 
survival distance of the nucleons is defined as 



m = 



&z 
H(zj 



H(z)(l + z) 



7i — 4 



(13) 



For ji ~ 2 and for those evolution scenarios T-L that we 
have considered so far in this paper, the effective survival 
distances range from 0.48/i?o ( n o evolution) to 2A/Hq 
(SFR). This agrees well with the relative ratio ~ 5 of the 
energy densities associated with lower neutrino limits in 
the proton-dominated scenario shown in Fig. [T] 



The relation ( 12 ) shows that as long as the maximal 



energy per nucleon is much larger than the pion produc- 
tion threshold in the CMB (i.e. i? max ^> AEqzk) and the 
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injection index is 7$ ~ 2 the main difference in the energy 
density of GZK neutrinos comes from the underlying evo- 
lution model, not by the inclusion of heavy elements. In 
principle, this factor can be large even for heavy nuclei if 
the sources have a strong evolution. The fact that typical 
CR models including heavy nuclei produce significantly 
less GZK neutrinos can be traced back to a low maxi- 
mal energy per nucleon and/or a weak evolution of CR 
sources. Note that the latter is an important ingredient 
of proton-dominated low-crossover models |28j . whereas 
CR models of heavy nuclei including more model degrees 
of freedom are less predictable w.r.t. the source evolution. 

Note that, ultimately, the inferred energy density w 7 
of the extragalactic diffuse 7-ray background in the 
GeV-TeV region constitutes an upper limit for the to- 
tal electro-magnetic energy from pion-production of UHE 
CR nuclei, see e.g. [24. An upper limit is given via the 
relation 



/ 1 4 



W GZK 



(14) 



Recent result from Fermi-LAT [26] translates into an en- 
ergy density of w 7 ~ 6 x 10~ 7 eV/cm 3 J2HH7]. Assuming 
an E~ 2 neutrino spectrum between energies E- and E + 
a numerical simulation gives a cascade limit of 



3 x 1CT 7 
log 10 (£+/£_) 



GeV cm 2 s 1 sr 1 



(15) 



This is only slightly lower than the estimate ( 14 ) for 
K n = 1. Cosmogenic neutrino fluxes that saturate this 
bound in the EeV region are already ruled out by Ice- 
Cube upper limits [11] , 



VI. DISCUSSION 



We have discussed in this article a simple procedure to 
derive lower limits on the cosmogenic neutrino flux. The 
limits are based on the observed spectrum and composi- 
tion of UHE CRs and depend on the unknown evolution 
of sources. For the case of a proton-dominance in the 
UHE CR data we show that ARA-37 should identify the 
flux of cosmogenic neutrinos after 3 years of observation 
if UHE CR sources follow the star formation rate. For 
the less optimistic (and less realistic) case of no source 
evolution it would require 10 years of observation. 

In the case of heavy nucleus dominance of the CR 
flux cosmogenic neutrino predictions are less optimistic. 



We can derive a lower limit in this scenario by track- 
ing the leading nucleus back to its source. Since photo- 
disintegration conserves the energy per nucleon of the in- 
teraction we can base our analysis on the observed num- 
ber of nucleons in UHE CRs, which depends on the ob- 
served mass composition. 

The dominant contribution to the cosmogenic neutrino 
flux is expected from the proton content in the UHE CR 
spectrum. We show in Fig. [2] two cases where we de- 
crease the contribution of protons to 10% and 1% at 
100 EeV and assume source evolution with the star- 
formation rate. Even this less optimistic case is in reach 
of ARA-37 after 5 years of observation. 

The prediction of cosmogenic neutrinos is very sensi- 
tive to the maximal CR injection energy per nucleon. 
If this is significantly larger than the GZK cutoff, even 
UHE CR scenarios dominated by heavy nuclei can pro- 
duce large fluxes of cosmogenic neutrinos. For flat spec- 
tra that are sufficiently close to E~ 2 the energy density 
of these optimistic GZK neutrino predictions depends on 
the cosmic evolution of the sources. 

All cosmogenic neutrino fluxes shown in this analy- 
sis are normalized to Auger data. The spectra observed 
with HiRes and the Telescope are in general larger, which 
could be a result of an overall systematic energy shift by 
20 — 30%. This corresponds to an upward shift of up to a 
factor 2 of the energy density E^ R Jcr(^cr)- Hence the 
lower limits shown in Figs. [T] and [2] should be similarly 
scaled upward. 

Finally, we would like to stress that the present anal- 
ysis does not take into account statistical uncertainties 
of the CR data. However, the method can be easily ex- 
tended to this case. In Ref. [21] we have shown that an 
actual fit to HiRes data assuming a proton power-law in- 
jection in the sources is statistically consistent with cos- 
mogenic neutrino fluxes that exceed the minimal bound 
by up to an order of magnitude and are in reach of the 
IceCube detector. 
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